QCD phase diagram at finite baryon and isospin chemical potentials 
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The phase structure of two-flavor QCD is explored for thermal systems with finite baryon- and isospin- 
chemical potentials, /ib and /iiso, by using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model. 
The PNJL model with the scalar-type eight-quark interaction can reproduce lattice QCD data at not only /liso ~ 
^B = but also /iiso > and fiB ~ 0. In the /liso-fJ-B-T space, where T is temperature, the critical endpoint 
of the chiral phase transition in the fiB-T plane at piso=0 moves to the tricritical point of the pion-superfluidity 
phase transition in the fiiso-T plane at /iB=0 as /Xiso increases. The thermodynamics at small T is controlled by 
Vct^+vt^ defined by the chiral and pion condensates, a and n. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

The phase diagram of quantum chromodynamics (QCD) is 
the key to understanding not only natural phenomena such as 
compact stars and the early Universe but also laboratory ex- 
periments such as relativistic heavy-ion collisions. Quantita- 
tive calculations of the phase diagram from the first-principle 
lattice QCD (LQCD) have the well-known sign problem when 
the baryon chemical potential (/ie) is real [1]; here, ^b is re- 
lated to the quark-number chemical potential /iq as /j,b = 3/Xq. 
Several approaches have been proposed so far to circumvent 
the difficulty; for example, the reweighting method ^, the 
Taylor expansion method |3|] and the analytic continuation 
from imaginary fiq to real /iq |3-Hl- However, those are still 
far from perfection particularly at fiq/T > 1, where T is tem- 
perature. 

As an approach complementary to LQCD, we can con- 
sider effective models such as the Nambu-Jona-Lasinio (NJL) 
model yiHlJl] and the Polyakov-loop extended Nambu- 
Jona-Lasinio (PNJL) model [|13M33|1 . The NJL model de- 
scribes the chiral symmetry breaking, but not the confine- 
ment mechanism. The PNJL model is extended so as to treat 
both the mechanisms 1 14] approximately by considering the 
Polyakov-loop in addition to the chiral condensate as ingredi- 
ents of the model. 

In the NJL-type models, the input parameters are usually 
determined from the pion mass and the pion decay constant 
at vacuum (/iq — and T — 0). Some of the models have 
the scalar-type eight-quark interaction. The strength of the in- 
teraction is adjusted to LQCD data at finite T ll3lll . since the 
sigma-meson mass at vacuum related to the interaction has 
a large error bar and then ambiguous |34]. It is then highly 
nontrivial whether the models predict properly dynamics of 
QCD at finite /iq. This should be tested from QCD. Fortu- 
nately, this is possible at imaginary /iq, since LQCD has no 
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sign problem there. In Ref. |]32]1 . it was shown that the PNJL 
model can reproduce LQCD data at imaginary fj,q. QCD has 
the Roberge-Weiss periodicity and the Roberge-Weiss transi- 
tion 1351] in the imaginary /iq region, because of the extended 
Z3 symmetry I13ll - [33ll . The PNJL model can reproduce these, 
since it has the symmetry |3lt]. In the real /iq region, as an 
important result, a phase diagram is predicted by the PNJL 
model with the parameter set |32] determined from the LQCD 
data at imaginary /iq. The PNJL prediction shows that in the 
/iq-T plane at /ii=0 there appears a critical endpoint (CEP), 
that is, a second-order critical point where a first-order chiral 
phase-transition line terminates. 

A similar test of the PNJL model is possible for finite 
isospin-chemical potential (//iso) [36], since LQCD has no 
sign problem there; for later convenience, we use the "mod- 
ified" isospin-chemical potential /ii — /iiso/2 instead of 
/iiso- LQCD data are available for both real JSTT and imag- 
inary |38, 39] /iiso. The PNJL model has already been applied 
to the real fj,i 1231 l24ll and the imaginary /ij case [40] with 
success in reproducing the LQCD data. The PNJL calculation 
at real iii shows that in the /ij-T plane at /iq = there exists 
a first-order pion-superfluidity phase-transition line connected 
to a second-order pion-superfluidity phase-transition line; the 
connecting point is a tricritical point (TCP) by definition. Real 
/ii dependence of QCD phase diagram is investigated in other 
models such as chiral perturbation theory 136] , the strong cou- 
pHng QCD nil] and so on IH. 

The CEP in the /iq-T plane at /ij = is important as a good 
indicator of the chiral and deconfinement phase transitions in 
collider experiments at Helmholtzzentrum fiir Schwerionen- 
forschung GmbH (GSI), Super Proton Synchrotron(SPS) 1,41 
H, Relativistic Heavy Ion Collider (RHIC) lH HI] and 
LHC ll47tl . In the measurements, /ij is not zero generally. It 
is then interesting to see how critical points such as CEP and 
TCP are located in the /iq-/ii-T space. 

In this paper, we draw the phase diagram of two-flavor 
QCD in the fii-^^-T space by using the PNJL model. Fol- 
lowing our previous paper ll32[l . we introduce the scalar- type 
eight-quark interaction to reproduce LQCD data on thermal 
systems at iiq — jii = 0. The scalar-type eight-quark interac- 
tion is a next-to-leading order coiTection in the power count- 
ing rule based on mass dimension. First, we will show that 



the PNJL model with the parameter set thus determined also 
reproduces LQCD data on thermal systems at /ii > and 
/iq = 0. After confirming the reliability of the present PNJL 
model, we will predict locations of CEP and TCP in the /ij- 
/iq-T space. 

In Sec. nil the PNJL model is recapitulated. In SecHUl the 
PNJL calculation is compared with LQCD data for thermal 
systems at /ij > and /iq ~ 0, and the phase diagram is ex- 
plored in the fii-^q-T space. Properties of the susceptibilities 
near CEP and TCP are analyzed. Section |IV] is devoted to a 
summary. 



II. PNJL MODEL 



The two-flavor PNJL Lagrangian in Euclidean space-time 



C = qhyD" - "fijl + mo)q + Gs [(qq)'^ + (qijbTq)'^] 
+ Gsiiqqf + {qil5fqf]' -KmA],mr ,T), (1) 

where D" ^ d" + iA" and A" = S^gA^^^ with the gauge 
field A^, the Gell-Mann matrix Xa and the gauge coupling 
g. In the NJL sector, Gs denotes the coupling constant of 
the scalar-type four-quark interaction and Gs stands for that 
of the scalar-type eight-quark interaction 11321 l48l 14911 . The 
Polyakov-potential U, defined in (fTsT l. is a function of the 
Polyakov-loop ^ and its Hermitian conjugate <!>*. 

The chemical potential matrix fi is defined by /i = 
diag(/itj , lid) with the w-quark (d-quark) number chemical po- 
tential flu ilJ'd), while rho = diag(mo, mo). This is equivalent 
to introducing the baryon and isospin-chemical potentials, /ie 
and /iiso, coupled, respectively, to the baryon charge B and to 
the isospin charge I3: 



A = MqTo + I^ITS 



(2) 



with 



/i„ + fid Mb fJ-u — Mrf Miso ,„, 



where tq is the unit matrix and t^ (i = 1, 2, 3) are the Pauli 
matrices in flavor space. Note that fiq is the quark chemical 
potential and fij is half the isospin-chemical potential (fiiso)- 
In the limit of toq = fii — 0, the PNJL Lagrangian has the 
SUl{2) X S'C/r(2) X [/v(l) X SUciS) symmetry. For mo ^ 
and/ii 7^ 0, itisreducedtoC/i3(l) x C/v(l) x SUd^). 

The Polyakov-loop operator 3 and its Hermitian conjugate 
^t are defined as 



S = —TrL, <pt = —TrL\ 

N ' TV ' 



with 



L(x) — Pexp 



I I dT^4(x, r) 




(4) 



(5) 



where V is the path ordering and A4 = iA'^. In the PNJL 
model, the vacuum expectation values, <P — {<P) and <P* = 



{'P''), are treated as classical variables. In the Polyakov gauge, 
L can be written in a diagonal form in color space ll§\: 

L = e^'3('^3A3+</.8As) ^ diag(e*'''^%e*'''^%e''^'^'=), (6) 

where (j)a ^ (j)3 + 08/\/3, (/ifc = -(/13 + ^s/V^ and (/i-c = 
-(0a + 06) = -208/ V3. 

The Polyakov-loop <P is an exact order parameter of the 
spontaneous Z3 symmetry breaking in the pure gauge theory. 
Although the Z3 symmetry is not exact in the system with 
dynamical quarks, it still seems to be a good indicator of the 
deconfinement phase transition. Therefore, we use <P to define 
the deconfinement phase transition. 

The spontaneous breakings of the chiral and the 1/13(1) 
symmetry are described by the chiral condensate a = (qq) 
and the charged pion condensate ll23ll 



T2'- 



^±iip 



{qij5T±q), 



(7) 



where t± = (ti ± 1x2) /\/2. Since the phase ip represents the 
direction of the [/i., (1) symmetry breaking, we take ip — for 
convenience. The pion condensate is then expressed by 



vr = {qijr-,Tiq). 



(8) 



Making the mean field (MF) approximation 111 iL 12311 . one can 
obtain the MF Lagrangian as 



C 



MF 



qijuD" - 74/i + Mtq + Ni^5Ti)q 
~Gs[a^ + TT^] - 3Gs{a^ 



7:^)^-U 



(9) 



with 



M 



mo - 2Gsa - AGsa{a^ 



N ^ -2Gs7r ~ 4G87r(cr2 + TT^). 
Performing the path integral in the PNJL partition function 



(10) 
(11) 



JPNJL 



DqDq exp 



d-'xC 



MF 



(12) 



we can get the thermodynamic potential i7 (per unit volume). 



[2 



-rin(ZpNjL)/T^ = -2^ /" 



(27r)3 



3Mp) 



+ ^ In [1 + 3(<? + <^*e-^^"^(P))e-^^.^(P) + g-^/^^rCp)] 



1 

G,[(t2 



In [1 + 3(<Z'* + <I>e-^''tip))e-f^Etip) + ^-^PE+ip)^ 
^]+3Gs{(J^+TrY+U (13) 



with i?±(p) = E±{p) ± /iq, where 

S±(p) = ViE{p)±f^i)' + N^ 



(14) 



for E{p) — y^j^ + AP. On the right-hand side of (fT3] i. only 
the first term diverges, and it is then regularized by the three- 
dimensional momentum cutoff A 11141 USll . 



We use U of Ref. UJ] that is fitted to LQCD data in the 
pure gauge theory at finite T IdOIdIII : 



U = T^ 



aiT) 



<P*^ 



h{T) ln(l - 6'^<?* + 4(<?-^ + 'P*") - 3{$$*y) 



(15) 



'To\3 



a{T) = ao + a,{^)+a,{^), 6(T) = 63(f) 

(16) 

where parameters are summarized in Table |l] The Polyakov 
potential yields a first-order deconfinement phase transition at 
T = To in the pure gauge theory. The original value of Tq 
is 270 MeV determined from the pure gauge LQCD data, but 
the PNJL model with this value of To yields a larger value of 
the pseudocritical temperature Tc at zero chemical potential 
than Tc = 173 ± 8 MeV that the full LQCD simulation fH- 
|54l predicts. Therefore, we reset Tq to 212 MeV [32] so as to 
reproduce the LQCD result. 



ao 



ai 



0-2 



3.51 



-2.47 15.2 



-1.75 



TABLE I: Summary of the parameter set in the Polyakov-potential 
sector determined in Ref. 1191 . All parameters are dimensionless. 



The classical variables X — <P,<P* 
by the stationary conditions 

df2/dX = 0. 



a, and tt are determined 



(17) 



The solutions to the stationary conditions do not give the 
global minimum of il necessarily. There is a possibility that 
they yield a local minimum or even a maximum. We then have 
checked that the solutions yield the global minimum when the 
solutions X{T, /iq, /^i) are inserted into ( fTSI ). 

In this work, first-order transitions are defined by (approx- 
imate) order parameters, c, vr and <1> in their discontinuities. 
When the susceptibility of one of the order parameters di- 
verges, we regard it as a second-order transition of the order 
parameter. For crossover, the pseudocritical point is deter- 
mined by a peak of the susceptibility. When the susceptibility 
has two peaks and it is not clear which peak should be taken, 
we do not plot a phase boundary to avoid the confusion. 

Table |ll] shows parameters in the NJL sector used in the 
present analyses. As shown in Ref. 1331, set A can reproduce 
not only the pion decay constant /tt = 93.3 MeV and the 
pion mass A/jr = 138 MeV at vacuum (T = fiq = ^i = 0) 
but also Tc = 173 ± 8 MeV JUSl at finite temperature 
(T > and ji^^ = Hi — 0). For this reason, we take this 
parameter set in this paper For comparison, we also use set 
B with no scalar-type eight-quark interaction. This param- 
eter set also reproduces the pion mass and the pion decay 
constant correctly, but not LQCD data at finite temperature 
(T > and /iq = /^i = 0). The sigma-meson mass M^ is 526 
(680) [MeV] for set A (B). 



Set 



Gs 



Gg 



mo 



A 



A 4.673 [GeV~2] 
B 5.498 [GeV-2] 



452.12 [GeV^*] 




5.5 [MeV] 631.5 [MeV] 
5.5 [MeV] 631.5 [MeV] 



TABLE II: Summary of parameters in the NIL sector. Here, To 
212 MeV for both the sets. 



III. NUMERICAL RESULTS 

The phase structure in the /ii-/iq-T space is explored by the 
PNJL model with the eight-quark interaction. 

A. Phase structure in the /xi-T plane at /iq = 

LQCD data are available in the ^i-T plane at /iq = 113711 . 
since LQCD has no sign problem there. In QCD, it is 
known |36] that at zero T, a second-order phase transition 
occurs at /ii — M.^/2 from the normal (tt — 0) to the pion- 
superfluidity phase (tt 7^ 0); this will be understood in sub- 
section IIIIBI also by using the PNJL model with the eight- 
quark interaction. The critical chemical potential /ic of the 
pion-superfluidity phase transition is /Xc = 0.57/a in LQCD 
calculation with a lattice spacing a, while it is /ic = Mt^/2 — 
69 [MeV] in the PNJL calculation. In the LQCD data, ^i is 
then normalized as /ic = 69 [MeV]. This makes it possible to 
compare the PNJL calculation with the LQCD data. 

First, we consider the normal phase by taking a case of 
/ii = 0.96/ic = 66 [MeV]. Figure [T] presents a and <? as a 
function of T /T^ where a is normalized by the value cto at 
zero T. LQCD data are plotted by plus (+) symbols with 
10 % error bar; LQCD data of Refs. [37] have only small er- 
rors, but we have added 10 % error that comes from LQCD 
data 1I53I1 on the pseudocritical temperature Tc at zero quark 
and isospin-chemical potentials. The PNJL result with the 
scalar-type eight-quark interaction (the thick-solid curve) is 
consistent with the LQCD data; note that the present model 
has no free parameter. If the scalar-type eight-quark inter- 
action is switched off from the PNJL model, the result (the 
thin-solid curve) deviates sizably from the LQCD data par- 
ticularly on a. This indicates that the scalar-type eight-quark 
interaction is inevitable. 

We use the dimensionless susceptibiUty matrix lll5ll26ll 



X 



c- 



(18) 



defined by the dimensionless curvature matrix 



C 



■K^ 



^TTIT ^TVa C-K^ C, 

^cTTT ^aa C(T^ ^a^ 
C^TT Ct^cr C^<^ C^^ 



\ ^'Ptt ^^a ^^^ C<p^ / 



T-^a 



•P-n 



T^a 

T-i/2, 



ipa 



T Qcpip T ii. 



ai> 



•P-P 



,(19) 



\ ^ '^PTT ^ ^^<P<T ^ 'i(pip J i'lptp / 



~1 T" 




1.2 1.4 



e 0.4 





Fig. 1: (color online). T dependence of (a) chiral condensate a and 
(b) Polyakov-loop <? at ^i = 0.96/ic = 66 [MeV] and pq = 0. Here, 
a is normalized by the value ao at vacuum and T is also normalized 
by Tc = 173 [MeV]. The thick (thin) solid curves represent the PNJL 
results with (without) the scalar-type eight-quark interaction; Lattice 
data (+) are taken from Ref. 1 37]. The lattice data are plotted with 10 
% error bar, since lattice calculations have 10 % error in determining 



where Q^y ~ d'^fi/dxdy for x,y = a,n,^,3. The sus- 
ceptibilities thus defined are dimensionless. For simplicity, 
we take the following shorthand notation: Xa ~ Xaa ■, Xt — 

Figure |2] presents the chiral and the Polyakov-loop suscep- 
tibility, Xcr and x^, as a function of T. The PNJL model 
with the scalar-type eight-quark interaction (the thick-solid 
curve) gives a better agreement with the LQCD data than 
the PNJL model without the scalar-type eight-quark interac- 
tion (the thin-solid curve). The present analysis for finite /ij 
is parameter free. Therefore, the reasonable agreement be- 
tween the PNJL model with the eight-quark interaction and 
the LQCD data indicates that the PNJL model with the eight- 
quark interaction is reliable. 

Next, we consider the pion-superfluidity phase by taking a 
case of /ii — \A^c = 96 [MeV]. Figure [3] presents <? and 



Fig. 2: (color online). T dependence of (a) chiral and (b) Polyakov- 
loop susceptibility at ^i = 0.96/ic = 66 [MeV] and /iq = 0. See 
Fig.[T]for the definition of lines and LQCD data. Since the suscepti- 
bilities of LQCD are obtained in arbitrary units, the magnitudes are 
then rescaled to fit the corresponding thick-solid curves, respectively. 



TT as a function of T/T^, where vr is normalized by the value 
ttq at zero T. Again, the PNJL model with the scalar-type 
eight-quark interaction (the thick-solid curve) is consistent 
with the LQCD data compared with the PNJL model without 
the scalar-type eight-quark interaction (the thin-solid curve). 
The PNJL calculation on tt shows that the pion-superfluidity 
phase transition is of second order there. 

Thus, the PNJL model with the scalar-type eight-quark in- 
teraction is consistent with the LQCD data, indicating that the 
model is more reliable than the original PNJL model without 
the eight-quark interaction. Figure|4]shows the phase diagram 
in the jii-T plane at /iq = 0. Panels (a) and (b) present results 
of the PNJL calculations with and without the eight-quark in- 
teraction, respectively. The thick-solid curve shows a first- 
order pion-superfluidity phase transition, while the dashed 
line indicates a second-order pion-superfluidity phase transi- 
tion. A meeting point between the two lines is a tricritical 
point (TCP) by definition. The dot-dashed (dotted) line stands 
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Fig. 3: (color online). T dependence of (a) Polyakov-loop and (b) 
pion condensate at fii — 1.4/^c ~ 96 [MeV] and /iq = 0. Here, n is 
normalized by the value ttq at zero T. See Fig.[T]for the definition of 
lines and the LQCD data. 
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Fig. 4: (color online). Phase diagram in the /ii-T plane at ^q = 
for the case (a) with and (b) without the scalar-type eight-quark 
interaction. The thick-solid (dashed) curve represents a first-order 
(second-order) pion-superfluidity phase transition. The dot-dashed 
(dotted) line means a deconfinement (chiral) crossover transition. At 
j-ii > M^/2, x* has two peaks, so we do not plot any deconfinement 
crossover transition line there. Meanwhile, the chiral crossover tran- 
sition line (dotted line) terminates at TCP. LQCD on the chiral and 
deconfinement crossover transitions are represented by a plus (+) 
symbol, while LQCD on the second-order pion-superfluidity transi- 
tion and the deconfinement crossover transition are shown by a cross 
(x) symbol. See Fig.[T]for more information on LQCD data. 



for a deconfinement (chiral) crossover transition. In panel (a), 
the two crossover transitions almost agree with each other In 
LQCD, meanwhile, the agreement is perfect, as represented 
by a plus (+) symbol with 10 % error bar LQCD data on 
the pion-superfluidity transition is also shown by a cross (x) 
symbol with 10 % error bar Comparing the PNJL results 
with LQCD data, we can confirm that the PNJL model with 
the eight-quark interaction is more consistent with the LQCD 
data than that without the eight-quark interaction. The loca- 
tion of TCP is (Aii,T) = (0.32 [GeV], 0.169 [GeV]) for the 
PNJL model with the eight-quark interaction and (/ii,T) = 
(0.401 [GeV], 0.171 [GeV]) for the PNJL model without the 
eight-quark interaction. Thus, the eight-quark interaction is a 
sizable effect also on the location of TCP. 



B. Phase structure in the A*i-Mq plane at T = 

In the /ii-/iq plane at T = 0, the thermodynamic potential 
of the PNJL model is reduced to that of the NJL model: 






d^p 



(27r)3 

^2l 



3Gs{a^+nY. (20) 



When ^i < A/^/2, tt = and M « 330 MeV, so that E± = 
E±ni>AI-^^ 260 MeV. Hence, when^q < 260 MeV, 
n is reduced to 



n = GsO-2 + 3G8CT'' 



12 



_fp_ 

(27r)3 



E{p). 



(21) 



Therefore, J7 does not depend on /i^ and /ij for /ij < Mt^/2 
and /iq < 260 MeV, indicating that no phase transition occurs 
there. In other words, there is a possibiUty that a chiral phase 
transition takes place when ^^ > 260 MeV. This is reahzed, 
as shown later in Fig.|6] For /ij = Mt,/2, more careful discus- 
sion is necessary, since it is a boundary of the normal phase in 
which TT ~ Q. This is discussed below. 

In the normal-phase region at /Xq < 260 MeV, the curvature 
of Q in the 7r-direction is obtained by (l 



Over the normal and pion-superfluidity phases, we have 



9? 



with 



20-77 — 48G-7r 



(27r)3 £;(p)2 _ ^2 



^^-"2 9^- 



(23) 



Thus, /(/ii) does not depend on jiq. Here, an effect of the 
eight-quark interaction appears only through M and G-j^. 

In vacuum (T = /iq = /tj = 0), the RPA function with 
external momentum (go 7^ 0, q = 0) is illl[l2ll27ll 



2Gr7I- — 48Gr7r 



(2^)3 i?(p)2 _ ^2/4 



/(f), (24) 



and the pion mass M-^ is determined by the condition 
f{M^/2) = 0. We then find for fii = M^/2 that 



9^ 



= f^^U 



V 2 



0. 



(25) 



Thus, the curvature of i7 at /ij = M^/2 is zero in the tt direc- 
tions, indicating that a second-order pion-superfluidity phase 
transition takes place at /ij = M.,^/2 when /Xq < 260 MeV. 
This point will be confirmed later in the phase diagram of 
Fig. |6la) where the second-order pion-superfluidity phase 
transition line (dashed fine) is a straight line. 

Next we consider both regions of /ij < A/,r/2 and jii > 
AU/2. Figure IS shows |7r|, \a\ and R = ^/IF^TN^ as 
a function of /ij and /iq. The pion condensate tt is zero at 
/ii < i\f7r/2, but nonzero at /ii > Mt;/2, as expected. There- 
fore, the former region is the normal (/3-symmetric) phase 
and the latter region is the pion-superfluidity (/a -symmetry 
broken) phase. The order parameter |cr| of the chiral sym- 
metry is almost constant in the normal phase but goes down 
in the pion-superfluidity phase. The parameter R is almost 
constant over the two phases, when /iq < 200 MeV. When 
/iq > 200 MeV, R has a discontinuity in the /iq direction. 
Thus, the /Xq dependence of R at finite /ij is similar to that at 
/ii = over a wide range of /ij. 

In the limit of niQ = and /ii = 0, the chiral symmetry 
is an exact symmetry. In this situation, the thermodynamic 
potential of ( fT3] l is a function of R, and i? is a function of 
Vc^ + 7r2. This means that R or Vc^ + '^'^ is an order pa- 
rameter of the chiral symmetry. When toq and/or /ij is finite, 
the chiral symmetry is not an exact symmetry anymore. How- 
ever, the fact that the /iq dependence of R at finite /ij is similar 
to that at /ii = means that the chiral symmetry is preserved 
with good accuracy. This is understood as follows. 



E± = ^{E±fxi)^ + N^ 



= \P 



i?2±2£;/ii + /i2. 



(26) 



As shown in Fig. IH i? is about 330 MeV at /iq < 200 MeV 
and III < A — 631.5 MeV. In the region, E± is well approx- 
/ij , because p^ 



imated by ^/Jp~^RF~^Jj^ , because p^ + i?^ + /ij ^ 2E^i. 
When /iq > 200 MeV, R is small and hence E± is approx- 

/(/ii), (22) imated by y/p^ + fif . Therefore, i7 is a function of R or 
-\/o'2 + 7r2 with good accuracy; here, note that R sa y/a'^ + n'^ 
because of niQ <C R- Thus, when T is small, the thermody- 
namics at finite /iq and /ij is controlled by an approximate or- 
der parameter it! of the chiral symmetry over both the normal 
(tt — 0) and the pion-superfluidity (tt 7^ 0) phase; the chiral 
symmetry is spontaneously broken when R is finite, while it 
is restored when R is zero. When T > T^, R is not large any 
more. Hence, a and tt work independently there, as shown 
later in subsection llllDI 

Figure |6] presents the phase diagram in the /ii-/iq plane at 
T — Q. When T — Q, the system is in the confinement phase 
because <? = there. So we consider the chiral and pion- 
superfluidity transitions only here. On the solid line, the first- 
order chiral and pion-superfluidity transitions coexist, while 
on the dot-dashed line only the first-order chiral transition 
takes place. The dashed line represents the second-order pion- 
superfluidity transition. In panel (a) where the eight-quark 
interaction is taken into account, the solid, dot-dashed and 
dashed lines meet at a point. This is a TCP, because the pion- 
superfluidity transition changes the order from first order to 
second order there, while the chiral transition keeps first or- 
der. Thus, there is no CEP in the /ii-/iq plane at T = 0. In 
panel (b) where the eight-quark interaction is switched off, the 
endpoint of the dot-dashed line is a CEP and a meeting point 
of the solid and dashed lines is a TCP by definition. Compar- 
ing the two panels, we can see that the eight-quark interaction 
changes the phase diagram qualitatively. 

Recently, it was shown in Ref. ||55|l that the A dependence 
of J7 may change the order of the phase transition in the mean 
field level. We then investigate the A dependence of the phase 
diagram in the /ii-/iq plane at T = 0. In this procedure, we 
consider the four-quark and eight-quark interactions only. The 
parameters of the PNJL model are determined for each A so as 
to reproduce the pion decay constant /^ — 93.3 MeV and the 
pion mass M^, = 138 MeV at vacuum and Tc = 173 ± 8 MeV 
at finite temperature with no /iq and /ii 154- 15411 ; note that T^ 
is a much stronger constraint on Gs than M„, since M„ has a 
large error bar [32]. This parameter fitting is exactly the same 
as that in Sec.|II]to determine the parameter set A. 

We vary A from 573 to 651.5 MeV. The upper and the 
lower bound of A are determined as follows. The QCD sum 
rule yields the lower and the upper bound of \cf\ as |cr| = 
(225 ± 25MeV)3 JMIH. The absolute value of the chiral 
condensate, |cr|, increases as A goes up, and reaches the upper 
bound of \g\ when A = 651.5 MeV. Thus, A = 651.5 MeV 
is the upper bound of A. Meanwhile, the lower bound of 
A is determined by not the lower bound of |cr| but the fact 



that no parameter set can reproduce /tt = 93.3 MeV and 
M^ = 138 MeV simultaneously when A < 573 MeV. Al- 
though this fact is found by numerical calculations, it can be 
understood with reasonable approximations. At zero temper- 
ature, the thermodynamic potential of the PNJL model is re- 
duced to that of the NJL model, as shown in ( l20l i. In the NJL 
model, the pion mass is obtained by 



Set 



Gs 



G; 



8 



mo 



A 



Ml = 



{M -mo)MI{M,M^) 



(27) 



with 



I{M,M^) 



8NfNc 
27r2 



p^dp 



/o v^p2^_^2[4(p2 + ^/2) _ ^/2] ' 

(28) 
where Nf and Nc are the numbers of colors and flavors, re- 
spectively, and Nf = 2 and A^c = 3 in the present case. Since 
M > Too in (|27]l, and AP > M^ in ^^, we then neglect 
Too in (l27b and Af^r in (|28] | in order to understand the math- 
ematical structure of ( |27] | and ( |28] |. Using the approximate 
equations and the Gell-Mann-Oakes-Renner relation, we have 



NfN, 2 



In 



1 + vr 



vr 



J 7T 

yl2' 



= ^, (29) 



where x — Mj A. The left-hand side of ( |29] l has a max- 
imum at X = 0.97, while the right-hand side increases 
monotonously as A goes down. This means that there exists 
a lower bound of A that satisfies ( |29] l. The lower bound is 
A = 573 MeV, although a obtained there is within the con- 
straint |cr| = (225 ± 25MeV)3 from the QCD sum rule. 

Table Hill presents three parameter sets. A, A' and A", ob- 
tained by the above procedure. Set A is the original parameter 
set mentioned in Sec. [Ill set A' is an example of the parameter 
sets near the lower bound of A, and set A" is the parameter set 
at the upper bound of A. The value of A in set A" is slightly 
larger than that in set A. This indicates that set A" yields qual- 
itatively the same phase diagram as set A. Actually, we have 
confirmed this with numerical calculations. Meanwhile, the 
phase diagram calculated with set A' is shown in Fig.|2] The 
phase structure shows no qualitative difference from the re- 
sult of set A in Fig.|6ja), although the first-order chiral transi- 
tion line (dot-dashed line) and the pion-superfluidity phase- 
transition line (solid line) are slightly shifted down by de- 
creasing A. Furthermore, we have confirmed that the phase 
diagram does not change qualitatively near the lower bound. 
Thus, the order of the phase-transition is not changed by vary- 
ing A in the range 573 < yl < 651.5 MeV. As a property 
of the parameter sets near the lower bound of A, Gg is quite 
large. This means that the higher-order multiquark interac- 
tions than the eight-quark interaction may not be negligible 
there. However, this sort of analyses is beyond the scope of 
the present work. 



C. Phase structure in the yi^-T plane at /ii = 

The phase diagram in the /iq-T plane at /ij = is shown 
in Fig. [8] The solid curve shows a coexistence line of first- 



A' 5.755 [GeV^2] 1264.2 [GeV"*] 5.77 [MeV] 580 [MeV] 
A 4.673 [GeV-^] 452.12 [GeV"*] 5.5 [MeV] 631.5 [MeV] 
A" 4.295 [GeV-^] 351.32 [GeV^] 5.31 [MeV] 651.5 [MeV] 

TABLE III: Cutoff dependence of parameters. Here, To = 203 MeV 
for the set A, To = 212 MeV for the set A and To = 217 MeV for 
the set A". 



order chiral and deconfinement phase transitions that ends at 
(^q, T) = (178 [MeV], 152 [MeV]). This point is a CEP by 
definition and is known to be of second-order [8, 10]. In gen- 
eral, once a first-order phase transition takes place for some 
order parameter, the discontinuity propagates to other order 
parameters unless the parameters are zero IBSIdST . The coex- 
istence between the first-order chiral and deconfinement phase 
transitions shown in Fig.[8]is a typical case of the coexistence 
theorem. The dot-dashed (dotted) line stands for a crossover 
deconfinement (chiral) transition. The crossover chiral and 
deconfinement transitions almost coincide with each other and 
end at the CEP. Thus, a CEP exists in the present model. This 
CEP survives, even if the eight-quark interaction is switched 
off. In the case of no eight-quark interaction, the CEP in the 
/iq-T plane at /ij = moves to a CEP in the /ii-/iq plane at 
T = of Fig. |6jb), as /ij increases from zero. This behav- 
ior of CEP is changed a lot by the eight-quark interaction, as 
shown later in Fig.fTOl 

In principle, the Polyakov-potential U depends on /iq as a 
consequence of the backreaction of the Fermion sector to the 
gluon sector. Particularly, the /iq dependence of the parameter 
To in U is important and estimated by using renormalization 
group arguments ll22ll : 



7o(Mc 



) = Ti-e "oM^q) 



(30) 



for 6(/iq) = 29/(67r) - Z^^iln-KT^) with ao = 0.304 and 



Tr = 1.770[GcV]. Figure |9] shows effects of To(/iq) on 
the phase diagram in the /iq-T plane at /ii = 0. Compar- 
ing this figure with Fig. [8] we can see that the effect dose 
not yield any qualitative change, but the location of CEP is 
moved from (fi^, T) = (178 [MeV], 152 [MeV]) to (/iq, T) = 
(187 [MeV], 130 [MeV]). At small T, the effect becomes neg- 
ligible, since U itself tends to zero as T decreases. 



D. Phase structure in the /ii-^iq-T space 

Figure [TOlpresents the phase diagram in the ^i-ji^-T space. 
In this space, TCP and CEP emerge not at points but on lines; 
precisely speaking, CEP appears on lines CD and DA, while 
TCP does on lines GD and DA. Thus, CEP moves from point 
C to A via D as /ij increases from zero. Meanwhile, TCP 
moves from point A to G via D as /iq increases from zero. 

Line GE is a second-order pion-superfluidity transition line 
in the /ii-/iq plane at T = 0. A track of the line with re- 
spect to increasing T becomes an area GEAD. Hence, the 



0.03 



0.02 



0.01 







0.1 
0.2 Hq [GeV] 



|o| [GeV ' 



0.03 



0.02 



0.01 







0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




0.35 — ■ 

0.3 
0.25 

0.2 
0.15 

0.1 

0.05 





(a) 



1st order a 

1st order o and k 

2nd order k 

TCP 



0.1 



0.2 



0.3 



0.4 



0.5 



m [GeV] 



1 St order o 

1 St order CT and 7t 

2nd order k 

TCP 

CEP 



0.1 



0.2 



0.3 



0.4 



0.5 



m [GeV] 



R [GeV] — 




Fig. 5: (color online). Order parameters (a) n, (b) a and (c) i? as a 
function of ^i and /iq. 



pion-superfluidity transition is second order on the area. Sim- 
ilarly, a track of line FG (GH) with respect to increasing T 
becomes an area FGDC (GHBAD). In area FGDC, the chiral 
and deconfinement transitions are of first order, while the pion 
condensate is zero. In area GHBAD, all the chiral, deconfine- 
ment and pion-superfluidity transitions are of first order. The 
two areas smoothly connect to each other, indicating that the 
thermodynamics in these areas are controlled by R. Properties 
of lines and areas in Fig.[TO]are summarized in TablellVI while 
locations of points in Fig. [TO] are summarized in Table IV] 



Fig. 6: (color online). Phase diagram in the fii-fiq plane at T = 
for the case (a) with and (b) without the eight-quark interaction. 
The solid line represents a coexistence line of first-order chiral and 
pion-superfluidity phase transitions, while the dot-dashed line shows 
a first-order chiral phase-transition line. The dashed line stands for a 
second-order pion-superfluidity phase transition. 
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Fig. 7: (color online). Phase diagram in the Hi-Hq plane at T = 
calculated with set A'. See Fig.|6]for the definition of lines. 



1st order 

cross over o 

cross over <I> 

CEP « 



0.05 



0.1 0.15 

(1. [GeV] 




T[GeV] 0.1 




0.1 
2 ^^q I'^SV] 



Fig. 8: (color online). Phase diagram in the /iq-T plane at ^i = 0. 
The solid line is a coexistence line of first-order chiral and decon- 
finement phase transitions. The dashed line stands for the chiral 
crossover transition, while the dot-dashed line does for the decon- 
finement crossover transition. Here, the eight-quark interaction is 
taken into account in the PNJL model. 
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Fig. 9: (color online). Effect of /iq-dependent To on the phase dia- 
gram in the /iq-T plane at /ii — 0. In the PNJL calculation with the 
eight-quark interaction, parameter To is replaced by /.iq-dependent 
parameter To(/iq). See Fig.[8]for the definition of lines. 



Figure [TT] presents the chiral susceptibility Xcr, the 
Polyakov-loop susceptibility X4> ^nd the pion susceptibil- 
ity Xtt as a function of /iq for the case of (/ij , T) — 
(0.075 [GeV], 0.140 [GeV]); these are plotted by the solid, 
dashed and dotted curves, respectively. The /Xq dependence 
of these susceptibilities correspond to a line parallel to the fiq 
axis in Fig. [TO] The susceptibilities Xa- and x* have peaks 
at the same position /iq = 187 MeV, indicating that the chi- 
ral and deconfinement transitions are second order there. This 
position corresponds to a point on line CD in Fig. [TO] Mean- 
while, Xtt has a peak at fiq = 173 MeV. This second-order 
critical point of the pion-superfluidity transition corresponds 
to a point on area ADGE in Fig.[TO] As an interesting feature, 
Xcr is discontinuous at /iq = 173 MeV. This property will be 
analyzed in Sec. lIIIE] 

Figure [T2I shows Xa, X* and x-^ as a function of /iq for the 
caseof (/ii,T) = (0.100[GeV],0.169[GeV]). All the suscep- 



Fig. 10: (color online). Phase diagram in the fii-fiq-T space. Here, 
the eight-quark interaction is taken into account in the PNJL model. 
Properties of lines and areas are summarized in Table II Vl while lo- 
cations of points are summarized in Table IVl 
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n $ 


CF 1st 


TT = 1st 


CD CEP 


TT = CEP 
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TT = <|i = 


AD CEP 


TCP CEP 


DG 1st 


TCP 1st 


GE 


2nd <1> = Q 


EA 


2nd 


GH 1st 


1st -I* = 



TABLE IV: Properties of areas and lines in Fig. [To] The phrase 
"1st" ("2nd") means that the phase transition either in the area or 
on the line is first (second) order. Blank means that no significant 
transition takes place there. 



tibilities have peaks at the same position /iq = 51 MeV, indi- 
cating that chiral, deconfinement and pion-superfluidity tran- 
sitions of second order take place simultaneously there. This 
critical point coiTesponds to a point on line DA in Fig. [TO] 
This is a TCP for tt and a CEP for a and (p. As an interesting 
feature, each of x<y and x<j has a kink at /iq ~ 51 MeV. This 
property will be analyzed in Sec. IIIIEI 

Now, the phase diagram in the /ii-/iq-T space is understood 
more precisely by considering the /iq-T plane at four values 
of /ij: each belongs to any of four regions, (i) /ij < ni{G) — 
AU/2, (ii) /ii(G) < ^il< (iiiD), (iii) iii{D) < /ii < fii{A) 
and (iv) ni{A) < /ij, where ni{X) is a value of /ij at point X. 
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TABLE V: Locations of points in Fig. 1101 
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Fig. 11: (color online). Chiral, pion and Polyakov-loop susceptibil- 
ities as a function of ^q at {ni,T) = (0.075[GeV], 0.140[GeV]). 
Here, the eight-quark interaction is taken into account in the PNJL 
model. These are represented by the solid, dashed and dotted, respec- 
tively. See Ref. 1 1 IJ for the definition of the susceptibilities. The x<y 
and Xrr are multiplied by 10""^ and 10~^, respectively, but x* is not 
multiplied by any factor. 
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Fig. 12: (color online). Chiral, pion and Polyakov-loop susceptibil- 
ities as a function of ^q at {fii,T) = (0.100[GeV], 0.169[GeV]). 
Here, the eight-quark interaction is taken into account in the PNJL 
model. See Fig.[TT]for the definition of lines. x<y ^rid Xtt are mul- 
tiplied by 1/20 and 10^*, respectively, but x* is not multiplied by 
any factor. 



1st order ■ 

1 St order o and O 

2nd order 7t 

cross orver o ■ 

TCP 

CEP 



0.15 
H„ [GeV] 



The fiq-T phase diagram in region (i) is essentially equal to 
that in the /iq-T plane at /ii = 0, i.e.. Fig. [8] since tt is always 
zero there. 

The /iq-T phase diagram in region (ii) is a bit more com- 
plicated, as shown in Figure [13] where /i/ = 75 MeV is taken 
as an example. In Fig.[T3j the thick-solid line ending at TCP 
stands for a coexistence line of first-order chiral, deconfine- 
ment and pion-superfluidity transitions. This is a natural re- 
sult of the coexistence theorem of the first-order phase transi- 
tion L33j 58J. Meanwhile, on the thin-solid line between TCP 
and CEP, first-order chiral and deconfinement transitions co- 
exist, but any first-order pion-superfluidity transition does not 
take place, because n is zero above the dashed line starting 
from TCP that represents a second-order pion-superfluidity 
transition. 

The fiq-T phase diagram in region (iii) is simpler than 
that in region (ii). Fig. [14] presents the fiq-T plane at 
/i/ = 100 MeV belonging to region (iii). As shown by the 



Fig. 13: (color online). Phase diagram in the /iq-T plane at fii = 
75 MeV. Here, the eight-quark interaction is taken into account in 
the PNJL model. 



thick-solid line, all the first-order chiral, deconfinement and 
pion-superfluidity transitions occur simultaneously there. A 
second-order pion-superfluidity transition and a crossover chi- 
ral transition occur on the dashed line start from a point shown 
by triangle. This point is a TCP for tt and a CEP for a by defi- 
nition. The point corresponds to a point on line DA in Fig.[TOl 

The /iq-T phase diagram in region (iv) is simple and easily 
imaginable from Fig. [TO] In this region, only a coexistence 
line of first-order chiral, deconfinement and pion-superfluidity 
transitions exists. 
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Fig. 14: (color online). Phase diagram in the fiq-T plane at fii = 
100 MeV. Here, the eight-quark interaction is taken into account in 
the PNJL model. The thick-solid line represents a coexistence line 
of the first-order chiral, deconfinement and pion-superfluidity transi- 
tions. On the dashed line, a second-order pion-superfluidity transi- 
tion and a crossover chiral transition occur simultaneously. 



E. Properties of susceptibilities 

Properties of the susceptibilities near the second-order 
pion-superfluidity transition line, CEP and TCP are investi- 
gated. 

For simplicity, we take the following shorthand notation for 
the curvature matrix C of ( fT9] i: 



C 




(31) 



where A = (cn-o-jC^.^ , Cj^.^), A^ is the transverse of ^, and 
the matrix K is expressed by 



^aa C(j^ ^(7^ 

^ ^ I C<fCT Ctpcp C^^ 

C/pa C^(p ^(73 



(32) 



As shown in JTSl ), fi is an even function of tt. Noting that 
cxY {.X,Y = ■K,a,<P,§) are proportional to d'^Q/dXpY; 
therefore, we can find that Cttt and c^tj for x^y = a,(l>,<P are 
TT-even, while Cxtt and c-^y for x,y = a,(l>,'P are 7r-odd. 

First, we consider the normal (tt = 0) phase including the 
second-orderpion-superfluidity transition line. Since tt = Oin 
this phase, the 7r-odd quantities Cx-k and c^y for x,y = a,<P,<P 
are zero. Therefore, C is reduced to 



C 



K 



(33) 



Equation ( l33T l shows the following properties. 



On the second-order pion-superfluidity transition line, 
the curvature c^^ in the tt direction is zero by definition 
of the second-order transition. Therefore, det[C] = 0. 
This indicates that Xtt = det[X]/det[C] diverges on 
the transition line, since det[i4r] is not zero in general. 



2. If a CEP of the chiral phase transition appears in the 
normal phase, the determinant dct[/ir] is zero at the 
CEP; see Ref. lUOtl for the details of this proof. Hence, 
Xa is divergent at the CEP because of det[C] = 
c^^det[A'] = 0. 

Properties 1 and 2 are understood clearly with numerical re- 
sults shown in Fig. [TT] The peak of Xtt at /iq = /i^ = 
173 MeV shows a second-orderpion-superfluidity transition, 
while the peak of Xa at /iq = /iq = 187 MeV does a CEP in 
the normal phase. Hence, the thermal system is in the normal 
phase (tt = 0) for /iq > /iq and in the broken phase (tt y^ 0) 
for /iq < /iq. Figure [15] shows dct[C] and dct[/ir] as a func- 
tion of /iq at (/ii,r) = (0.075[GeV],0.140[GeV]). It is found 
from this figure that det[C] ~ and det[/4r] 7^ at /iq = //q, 
while det[C] — det[K] = at /iq = /ijl. Thus, properties 1 
and 2 are confirmed to be true by the numerical results. 

Next, we consider the broken (tt 7^ 0) phase. Fig. [TTlis 
a good example. At /iq slightly smaller than /iq, tt is small, 
because tt — at /iq — /iq. Hence, any quantity can be 
expanded into a power series of tt. After the expansion, the 
TT-even quantities Cxy (x, y = (Tj^P, 3) are of order (7r)°, while 
the TT-odd quantities c^j, and Cxn (x, y = a,<P, 3) are of order 
(tt)^. The entry c^jr is of order (tt)^, as shown below. The 
stationary condition dTTI i for tt is rewritten into 







Ott 



dQ dn^ 
dn^ dn 



dQ 



2TT, 



and hence 



df2 



= 



(34) 



(35) 



because of tt 7^ 0. Expanding the 7r-even function J7 into a 
power series of tt^. 



[2 — y a„7r 

n 

one can see from ( 1351 ) that 

ai =0. 



(36) 



(37) 



Hence, Ctttt — T^d^fi/diidiT is of order (tt)^. Therefore, the 
matrix C is the following property in the broken phase: 

3. At /iq slightly smaller than /Xq, Ctttt is of order (tt)^, A 
and AJ- are of order (tt)^, and K is of order (7r)°. 

Now, we consider the reason why Xcr is discontinuous at 
/iq — /iq in Fig.[TT] The susceptibility x^ is expressed by 



Xrr = 



det[q ' 



(38) 



where A^cr is the cofactor of entry Cca in the matrix C. Prop- 
erty 3 indicates that both A^a and det[C] are of order (vr)^ in 
the broken phase at /iq < /iq, so that the left-hand limit of Xa 
as /iq approaches /iq is finite. As an important point, the tt- 
odd quantities contribute to this left-hand limit. Meanwhile, 
the TT-odd quantities are zero in the normal-phase at /iq > /iq. 
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so that they do not contribute to the right-hand limit of Aa-rr 
and det[C] as /iq approaches /i^. Thus, the right-hand limit 
of Xa is different from the left-hand limit of x„ . 

In Fig. [T2I all the susceptibilities, Xcr> X* ^"d Xtt, have 
peaks at the same position /iq = /i^ = 51 MeV. The di- 
vergence of x<T means that in dSST l. the denominator det[C] 
tends to zero faster than the numerator Aca as /iq approaches 
/iq from the left-hand side. There is no guarantee that such 
a strong damping of det[C] also happens in the right-hand 
limit, because 7r-odd quantities Cay and Cxtt are zero there. 
Actually, such a fast damping in the right-hand limit does not 
occur here, as shown by the numerical calculation. As /iq ap- 
proaches /iq, therefore, Xa- is divergent in the left-hand limit, 
but finite in the right-hand limit. This is the reason why Xa 
has a kink at /iq = 51 MeV. 

The fast damping of det [C] in both the left- and the right- 
hand limit happens only on point D in Fig. [TOl as shown be- 
low. Figure [T6] presents Xo-, X* and Xtt as a function of /iq 
at /ii = 0.08425 GeV and T = 0.136 GeV. All the suscepti- 
bilities diverge at /iq — 0.190 GeV. This peak corresponds to 
point D in Fig. [10] In this case, obviously, the susceptibilities 
have no kink. Therefore, det[C] tends to zero faster than Aa-a 
in both the right- and the left-hand limit. Point D is a meet- 
ing point of CEP and TCP. There is no guarantee that such a 
special critical point always happens. Actually, such a point 
does not appear if the eight-quark interaction is switched off, 
as shown in Fig.[T7l here, line CD (AG) represents CEP (TCP) 
of the chiral (pion-superfluidity) phase transition and there is 
no meeting point between CEP and TCP. 
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Fig. 16: (color online). Chiral, pion and Polyakov-loop susceptibili- 
ties as a function of /iq at (/ii, T) = (0.08425[GeV], 0.136[GeV]). 
Here, the eight-quark interaction is taken into account in the PNJL 
model. See Fig. [TT] for the definition of lines. The Xa ^rid Xtt are 
multiplied by 5 x 10^* and 2 x 10^^, respectively, but x* is not 
multiplied by any factor. 
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Fig. 17: (color online). Phase diagram in the /ii-/iq-T space calcu- 
lated with the PNJL model with the four-quark interaction only. See 
Fig.[TO]for the definition of lines and areas, except line CD represents 
CEP and line AG stands for TCP 



Fig. 15: (color online), /iq dependence of det[C] and det[_K'] at 
T = 0.14 [GeV] and m = 0.075 [GeV]. The solid (dashed) line 
stands for det[C] (det[i(']). Here, the eight-quark interaction is taken 
into account in the PNJL model. The det[C] is multiplied by 6 x 10^ . 



IV. SUMIVIARY 

Critical points such as CEP and TCP are important as in- 
dicators of the chiral, deconfinement and pion-superfluidity 
phase transitions in measurements at GSl, SPS, RHIC and 
LHC. In the measurements, /ij is not zero generally. We have 
then predicted the phase diagram of two-flavor QCD in the 



/ii-/iq-r space by using the PNJL model with the scalar-type 
eight-quark interaction. The PNJL model with the scalar-type 
eight-quark interaction is consistent with the LQCD data IBTIl 
in the fii-T plane at /tq = 0, while the original PNJL model 
without the scalar-type eight-quark interaction is not. 

In the /tq-/ii-T space, as shown in Fig.[T0] a CEP in the /Xq- 
T plane at fii — moves to a TCP in the fii-T plane /iq = 
as /ti increases. Meanwhile, the TCP in the /ii-T plane at 
/iq = moves to a TCP in the /iq-/ii plane at T = 0. When 
/ii < M.^/2, the pion condensate tt is zero and hence a CEP 
exists but any TCP does not. When M^/2 < /ij < 80 MeV, 
a CEP and a TCP exist separately. And when /ij > 80 MeV, 
they coexist. If the eight-quark interaction is switched off, a 
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CEP in the i-iq-T plane at ^i = moves to a CEP in the /iq-/ii 
plane atT — Q as fj,i increases; see Fig. [IT] Thus, the eight- 
quark interaction changes the QCD diagram qualitatively in 
the jjLq-fjLi-T space. 

When T is small, the thermodynamics at finite ^i and ^q 
is controlled by \/a'^ + tt^. The quantity Vo^ + ti^ is an ap- 
proximate order parameter of the chiral symmetry over the ly,- 
symmetric (tt — 0) and Is -symmetry broken (tt ^ 0) phases. 
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